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1. INTRODUCTION 

Construction of road embankment on the volt ground often leads to problems in the design stages, 
Soli sil is characterized as weak soil due to its igh compressibility and low shear strength, For ths reson, 
the soft ground should be weated to enhance the embankment stability. Various soft ground improvement 
methods cant be used, sich as the installation of piles, prefabricated vertical drains (PVDs) and 
stone columns. Installation of PVDs is a popular method at present duc to its effectiveness to accelerate 
‘consolidation and time-saving construction. 

‘Slope stability is the main factor contributing to the workability of road embankment (1) 
‘Traditionally, the method of limit equilibrium [2} and finite element [3] ure used to estimate slope stability 
Recent trends, artificial intelligence approaches have been well received by researchers in various 
fields of study. The application of the artificial intelligence approach is crucial in solving complex tasks to 
save time, energy and improve service quality [4]. It also has great potential in automating repetitive 
‘operations and work. 
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‘The artificial intelligence technology is one of the branches of science in making objects as 
intelligent as humans [5]. The important pat of artificial intelligence is machine leaning, which isthe ability 
ff the machine to lear to he smart. Conceptually, the machines that human will aced to learn must be 
‘eguipped witha lot of information or data to make them smarter than other machines [6 In other words, this 
‘machine will be equipped with thousands of altempls to complete a task, During this process, the machine 
‘will understand and learn concepts to solve the ask, 

“There are several approaches to artificial inlelligence approaches such as artificial neural networks 
(ANN), adaptive neuro-fuzzy inference system (ANFIS) und support vector machine (SVM). However, ANN 
‘s often used in engineering applications. ANN is a network designed to resemble the human brain intended 
to perform a specific task [7]. One ofthe important features of ANN model analysis is the number of neurons 
in the hidden layer. Inconsistencies in the number of neurons in the hidden layer cause the network to 
respond that diverges tothe target [SI 

‘A key aspect for producing an ANN prediction model with relatively high aceuracy is optimized 
based on the number of neurons [9], There ate various numbers of neurons reported from the results of 
previo studies. Due to the number of neurons inthe hidden layer as been reported to be inconsistent, there 
4 increasing concem among researchers [10-12], The issue of the number of neurons has been a 
‘controversial subject and widely debated in the field of arficial intelligence. Most previous studies have 
‘only focus on optimizing the earning algorithm approsch. Until recently, there was no reliable evidence that 
the number of neurons was determined to be used inthe ANN prediction model fr all engineering eases 

‘This paper secks to present an appropriate ANN model for predicting slope stability with high 
accuracy. The optimum number of neurons in the hidden layer contributing to the prediction aceuracy was 
{determined based om the trial and error method. In order to inaprove numerical simulation results, validation 
fof the smear effect parameter was performed by comparing the ld data. This study provides an exciting 
‘opportunity to advance our knowledge in improving the prediction aecuracy of ANN models 


2, NUMERICAL SIMULATION PROCEDURE 
‘A road embankment built on the soft ground treated with PVDs, as shown in Figure 1 is used as a 
case study. The embankment is built to a height of 2m while the slope is 1V: LH. A geotexile is placed on 
the soft ground surface as reinforcement to enbance embankment stability. The sand cushion was built on the 
soil surface tothe height of 0.5m. The soft soil heneath the embankment consists of three layers, namely soft 
silty clay, sity clay, and clayey silt. The location of the groundwater level (GWL) is 2.3m from the soft 
‘ground surface. The PVDs were installed in soft ground to «depth of 4m witha distance between cach other 
is Im, Two settlement gauges (SGI and SG2) were installed in the middle and toe of embankment on the soft 
‘ground surface to obtain measured data, Field measured data is observed for 177 days and collected weekly. 


Figure 1. Cross-section of soil profile and instruments installation location 


A series of laboratory tests were performed to determine the physical and engineering properties of 
the soft soil. The material parameters used in the numerical simulation are presented in Table 1, Numerical 
simulations were performed with software Plaxis 2D (version 8) by a plane strain model. The medium-mesh 
density consisting of 15-node triangular clements was used for numerical method analysis. The smearing 
elfect was considered in the consolidation analysis of soft ground with PVDs. The equivalent horizontal 
permeability f,) proposed by Hird etal. [13] was used to develop the plane strain model. According to the 
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literature, the smear effect of permeability ratio of 34 and 5 are used to very bye values that produce 
accurate numerical simulation output. The material models of the embankment and soft ground are Mohr- 
Coulomb (MC) and Soft Soil (8S), respectively. In addition, PVDs is modelled with drains function while 
_geotesiles isa geogrid function, 


Table 1. The numerical parameters of the embankment and soft ground, 


Fico sor mf 0 0 0 


3. NUMERICAL RESULTS AND DISCUSSION 

Tn order to provide the ANN dataset, validation of the smear effet parameter was performed by 
comparing the results of the sumerical simulation and measured. Correlation accuracy performance was 
assessed with the coefficient of determination (R*) based on the soil setlement data. The setlement curves 
produced at SGI and SG 2 has been illustrated in Figure 2a and b for 177 days. On the whole, the simulation 
results with 9 = 3 have shown considerably higher settlement than the other series. This may be because the 
patter of permeability distribution in the soft soil is not consistent. Soft soi is characterized as high soil 
‘Water content, and when forced to flow out through the PVDs, it will ereate an imbalance in the flow rate 
‘This has been confirmed by several researchers who have investigated and discussed this isue [14-16] 
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Figure 2. The relationship of soil settlement and time, (a) SGI. (b) SG2 


‘Table 2 presents the R? of settlement from numerical simulation and measured. Prom the data in this 
table, it s clear that showing the accuracy of the numerical simulation model of settlement with 7=4 yields 
the highest R? values in SGI and SG2, respectively, 0.9858 and 0.9614. This result is in agreement with the 
findings of [17], Interestingly, the output generated from the numerical method of ground settlement is 
considered good duc to the values of R* exceeds 0.9. This indicates thatthe findings support the literature 
‘where the permeability ratio of smear effect with values of 3, 4 and 5 ean be used for the consolidation 
analysis for soft ground treated with PVDs. 
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‘Table 2. The R? values obtained for numerical simulation 


4. ARTIFICIAL NEURAL NETWORK 

‘Anifcial neural network is a system that can acquire, store, and apply knowledge gained from 
experience. It can be used to organize and categorize information in ways that have proven effective in 
solving complex problems, dificult to understand or require additional information when solving with 
traditional computational methods [18]. Basically, ANN is composed of several processing units. Neural 
networks are developed based on the relationships between units that will respond in parallel to each input 
signal provided . To develop a neural network, the number of processing units must be determined at the 
input and output levels as well asthe type of node that corresponds tothe form of network to be used, 

Generally, nodes have three layers: the input node, the hidden node, and the result node. Each node 
receives input signals via weighted and responds tothe line weighted as output [19]. On the other hand, every 
relationship has its weight value. The weighted will be continuously updated when the learning process is 
performed on the network. Besides, it will always be adjusted so that the input and output behavior of the 
rework is loser to the infrastructure that provides the input value [20]. ANNs ean be presented in a more 
‘complex form where they contain a number of neurons that connect in which the output signal of one neuron 
will be th input signal to one or more neurons (21) 

‘As inthe biological nervous system, the relationship between the elements determines the function 
fof the network. It can be trained to perform specific tasks by modifying the value of the connector or 
‘weighting between elements. Typically, the neural network wis trained or modified to produce a specific 
input or output target specific. The modified neural network based on a comparison between the output and 
the target value so that the output value equals the value ofthe target network, 


‘8. Optimum number determination of neurons in a hidden layer 

‘More recent attention has focused on the optimum umber of neurons in the hidden layer (22-24) 
‘The small number of neurons will eause the network is notable to learn the relevance of the data. On the 
other hand, excessive numbers of neurons in the hidden layer interfere with network capability and cause the 
network to overfiting. Determining the number of neurons in hidden layer is complicated to perform [25] 
To date, there is no mathematical method that can be used to obtain the optimum number of neurons. The 
‘number of neurons inthe hidden layer can be obtained afte tying Various network structures because there 
is still no theory tha can be used to determine neurons number required in the hidden layer 

Generally. the numberof neurons in the hidden layer is determined using a tial-and-error method, 
In this study, the optimum number of neurons in the hidden layer was determined by the Wial-and-error 
method. The number of neurons in the hidden layer was evaluated ranging from one to 10 neurons. The 
‘number of neurons in the hidden layer is best to give the lowest eror. Thus, the optimum number of neurons 
‘was determined based on two performance indices, namely, the values of rot mean square error (RMSE) and 
the correlation coefficient 


1b, ANN models development 

"MATLAB Neural Networks Toolbox version R2010b was used to analyze input and output data for 
ANN models. Table 3 shows the input and output parameters used for ANN prediction. The outpa parameter 
is used asthe target value to evaluate the performance ofthe prediction 


‘Table 3. The input and output parameters for predictions 


‘The ANN model structure used in this study is three inputs and one output. The network structure 
used isa multi-layer network consisting of input layers, hidden layers, and output layers. The network used is 
feed-forward with tan-sigmoid transfer function in the hidden layer and linear transfer function in the 
‘output layer. 
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In ANN model analysis, data sets are randomly and automaticaly divided into three stages, 
‘namely training. validation and testing. The input and output vectors were divided into 70% (53 samples) for 
training, 18% (11 samples) for validation and another 15% (11 samples) for network testing, The flow chaet 
{oe predicting slope stability in the stady using the ANN model is shown in Figure 3, Training on the network 
‘was performed using a Levenberg-Marquardt backpropagation algorithm. The Levenberg-Marquardt 
backpropagation learning algorithm is designed specifically for squared error functions Its more likely to be 
used because ofits faster implementation than other types of backpropagation, such asthe conjugate gradient 
and resilient 


Figure 3, Flowchart slope stability prediction using the ANN model 


Performance index 

"The prediction performance of the ANN model yas evaluated by comparing the predicted output 
values with the larget values using the coefficient of determination (J), and the root means square error 
(RMSE). If the value of f° is one, it means that the output value is equal to the target value and are 
considered a good model prediction. The smaller RMSE values are closer to the target values, and if its 
values are zero, it means no error and indicates the accuracy of the prediction model is good. The definitions 
‘ofthe statistieal parameters used in this tady to evaluate the ANN models ure presented i (1) and (2) 
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where yy is the mean values, yoy, is the measured values, yj is the predicted values, and N 
is the dita number. 

‘The performance index used is capable of ensuring thatthe proposed model can provide a consistent 
level of accuracy in all periods. The advantages of using this statistical index are as performance indicators 
foe the proposed model and ensuring thatthe results of errs when assessing model performance are small 
and acceptable errors. Besides, the use of these performance indices can guarantee a consistent level of 
terror thal potentially provides the same level of error as the model is eximined for invisible data during 
the test period, 


Tat Ami Tatell, Val, 9, No, 2 une 2000, 256243 


Tot Arif Inell ISSN: 2252-4938, a ow 


PREDICTED RESULTS AND DISCUSSION 

In order to assess the errors that arise during training, validating and testing in the ANN model, K® 
and RMSE were used. Trial and error method to determine the best number of neurons inthe ANN network 
has been performed. Figure 4a and b show the relationship between the number of neurons with RMSE and 
HE, respectively. What is interesting in this graph is thatthe change inthe number of neurons clearly affects 
the prediction performance. It can be seen that prediction performance changes significantly fom the number 
fof neurons 1 t0 10, However, the subsequent increase in the number of neurons from 11 10 20 caused the 
prediction performance to decline, where the RMSE and X° values were relatively consistent. The use of a 
large number of neurons inthe hidden layer can lead to averiting. 
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Figure 4. Relationship between statistical index and number of neurons, (a) RMSE, (b) 


Overfiting is a major issue when designing neural networks because the actual data obtained is 
small relatively. If the numberof parameters used in a network i significantly larger relative to the data size 
used, the neural network will end to memorize rather than generalize the data properly. As a result, network 
performance during testing and applications will decline, The results show that the RMSE values are close to 
zero, and the R® valves are close to one for predicting factor of safety (FOS) found in neurons of 6, Thus, the 
‘number of neurons six was chosen for use in the ANN network to predict slope stability 

‘A comparative analysis between target and output was performed to predict the safety factor of the 
embankment slope. Figure 5 shows the relationship between target and output. It can be seen that the target 
‘and output data seem to overlap with each other. On the other hand, the resulting error between these dat is 
illustrated in Figure 6, where the error range is -0.05 to 0.0. These results suggest that large amount of data 
ig required while raining the network to produce good predictions. This is due to the neural network needs to 
lear the whole pattern, The total data used will alfect the performance of the network, especially in the 
testing and application stages 
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Figure 5. The relationship between FOS values calculated andl ANN 
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Figure 6. The err values between calculated and ANN 


6. CONCLUSION 

‘The main goal of the current study was to predict the slope stability for embankment built on the 
soft ground treated by PVDs. This study shows tht the neurons number in the hidden layer contributes to the 
prediction accuracy. These findings suggest that the optimum neuron number can be achieved by the trial and 
terror method for ANN training. The second major finding was that the permeability ratio ofthe smear effect 
is useful by applying a comparative analysis method between numerical simulation results and measured. 
‘The present study makes several noteworlhy contributions to the ANN model prediction for improving the 
accuracy and method of estimating the permeability ratio of smear effect in consolidation analysis. The major 
limitation ofthis stady is thatthe sample size used in ANN traning is small at 75. What is now needed is @ 
‘comparative study involving variations in the numberof samples to investigate its influence onthe prediction 
accuracy of ANN models, 
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